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I review the available empirical information on the equation of state of cold strongly 
interacting matter, as well as the prospects for obtaining new insights from the 
experimental study of gravitational waves emitted by neutron stars. 



1 Introduction 

The equation of state (EOS) is a nontrivial relationJinking the thermodynamic 
variables specifying the state of a physical system EJ. The best known example 
of EOS is Boyle's ideal gas law, stating that the pressure of a collection of 
N noninteracting, pointlike classical particles, enclosed in a volume fi, grows 
linearly with the temperature T and the average particle density n = N/Q. 

The ideal gas law provides a good description of very dilute systems. In 
general, the EOS can be written as an expansion of the pressure, P, in powers 
of the density (we use units such that Boltzmann's constant is Kb = 1): 

P = nT [l + nB{T) + n 2 C{T) + ...] . (1) 

The coefficients appearing in the above series, that goes under the name of 
virial expansion, are functions of temperature only. They describe the devia- 
tions from the ideal gas law and can be calculated in terms of the underlying 
elementary interaction. Therefore, the EOS carries a great deal of dynamical 
information, and its knowledge makes it possible to establish a link between 
measurable macroscopic quantities, such as pressure or temperature, and the 
forces acting between the constituents of the system at microscopic level. 

This point is best illustrated by the van der Waals EOS, which describes 
a collection of particles interacting through a potential featuring a strong 
repulsive core followed by a weaker attactive tail (see ng|]). At \U \/T « 1, 
Uo being the strength of the attractive part of the potential, the van der Waals 
EOS takes the simple form 

tlT 

P = - - an 2 , (2) 

1 — no 
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and the two quantities a and 6, taking into account interaction effects, can be 
directly related to the potential v(r) through 

r°° 16 
a = n \v(r)\ 2 r 2 dr , b — — nr^ , (3) 

J2r 3 

where 2r denotes the radius of the repulsive core (see fig.|l]). 
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Figure 1. Behavior of the potential describing the interactions between constituents of a 
van der Waals fluid (the interparticle distance r and v(r) are both given in arbitrary units). 




In spite of its simplicity, the van der Waals EOS describes most of the 
features of both the gas and liquid phases of the system, as well as the nature 
of the phase transition. 

The EOS of strongly interacting matter reflects the complexity of the 
underlying dynamics, that makes the phase diagram extremely rich. As an 
example, fig.^j shows tha-phase diagram of charge neutral strongly interacting 
matter in /3-equilibrmrrn. 

In this talk, I will review the available empirical information on the EOS of 
strongly interacting matter at T ~ 0, and the extent to which this information 
can be used to constrain theoretical models. The knowledge coming from 
nuclear systematics and neutron stars data is summarized in Sections 2 and 
3, respectively. Section 4 provides a short overview of the theoretical models 
of neutron star matter at supranuclear density, while Section 5 is devoted 
to the prospects for obtaining new insights from the experimental study of 
gravitational waves emitted by neutron stars. Finally, Section 6 summarizes 
the main results and states the conclusions. 
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Figure 2. Temperature vs baryon density phase diagram of charge neutral strongly interact- 
ing matter in /3-equilibrium. Hatched areas correspond to mixed phases of hadronic matter 
(HM) and quark matter (QM/QGP), as well as the nuclear liquid-gas (from ref.0)- 



2 Constraints on the EOS at T—0 from nuclear data 

Under standard terrestrial conditions, strongly interacting matter is observed 
in form of rather small, cold chunks: the atomic nuclei. 

Nuclei are self-bound systems consisting of Z protons and (A — Z) neu- 
trons. Their smallness is a consequence of the electrostatic repulsion between 
protons, limiting A to ~ 200, while the fact that they can be described as 
cold objects follows from the observation that thermal energies are negligible 
in comparison to the large proton and neutron Fermi energies. 

The body of data on nuclear masses provides a constraint on the density 
dependence of the energy per particle e = E/A at zero temperature, related 
to the EOS through 

^-'•-(iL-'dL- 

The A-dependence of the nuclear binding energy (i.e. the difference be- 
tween the measured nuclear mass and the sum of the constituent masses) is 
well described by the semiempirical formula 

S(Z, A) = a v A + a s A 2 / 3 + + a A (A - A a p -^ , (5) 
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whose terms account for nuclear interactions, taking place both in the interior 
of the nucleus and on its surface, electrostatic interactions between protons, 
isospin asymmetry and shell effects. 

The coefficient of the term linear in A yields the binding energy of sym- 
metric nuclear matter, an ideal uniform system consisting of equal number of 
protons and neutrons coupled by strong interactions only. The equilibrium 
density of such a system, no, can be inferred exploiting saturation of nuclear 
densities, i.e. the fact that the central density of atomic nuclei, measured by 
elastic electron-nucleus scattering, does not depend upon A for large A. 

The empirical equilibrium properties of nuclear matter are 

e = e(n = n , T = 0) = -16 MeV/A , n ~ .16 fm~ 3 . (6) 

In the vicinity of the equilibrium density e(n,T = 0) can be expanded accord- 
ing to (as we will only discuss the EOS at T = 0, the dependence upon T will 
be omitted hereafter) 

, \ 1 K (n- n Q ) 2 . . 

where 



n—no 



is the (in)compressibility module, that can be extracted from the measured 
excitation energies of nuclear vibrational states. Empirical values range from 
K ~ 200 MeV (corresponding to more compressible nuclear- matter, i.e. to a 
soft EOS) to K ~ 300 MeV (corresponding to a stiff EOS) 0. 

It has to be emphasized that the quadratic extrapolation of Eq. (ffi) can- 
not be expected to work far from equilibrium density. In fact, assuming a 
parabolic behavior of e(n) at large n (>> no) leads to predict a speed of 
sound in matter, v s , larger than the speed of light, i.e. 

(?)"5©> 1 - < 9 » 



regardless of the value of K. 

Eq.@ shows that causality requires 



dP 



< 1 , (10) 
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where e = E/Vt is the energy-density. For a noniteracting Fermi gas e oc ri 
implying (the equal sign corresponds to massless fermions) 



4/3 



' (7)^3 



(11) 



In presence of interactions the above limits can be easily exceeded. For ex- 
ample, modeling the repulsion between nucleons in terms of a rigid core leads 
to predict infinite pressure at finite density. 

The stiffest EOS compatible with causality is P = e, corresponding to 
(v s /c) = 1. Back in the early 60's Zel'dovich was able to show that the v s = c 
limit, corresponding to e oc n 2 , is indeed attained in a simple semirealistic 
theoryJn which nucleons are assumed to interact through exchange of a vector 
meson □. 

3 Constraints on the EOS at T—0 from neutron stars data 

At very large A ( ~ 10 57 ), gravity becomes strong enough to balance the re- 
pulsive interactions between nucleons, and nuclear matter can occupy macro- 
scopic regions of space. This is the situation occurring in the interior of 
neutron stars, compact astrophysical objects (typical values of mass and ra- 
dius are ~ 1.4 M Q and ~ 10 Km, respectively) whose central density largely 
exceeds n . The structure of a neutron star is schematically represented in 



Figure 3. Cross section of a neutron star. Note that the equilibrium density of nuclear 
matter corresponds to ~ 2.7 X 10 14 g/cm 3 . 



fig]3J. Note that the inner and outer crusts contain a comparatively small 



inner crust [~ 0.5 Km) outer crust (~ 0.3 Km) 

nuclei + n + e ^irzs^^^^z^s-^ nuclei + e 




pss2xl0 14 g/cm ; 



PS-4X10 11 g/cm 3 



^ R ~ 10 Km 
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amount of matter, most of the star mass being concentrated in the region of 
supranuclear densities (n > no). In addition, as all relevant Fermi energies are 
much larger than the typical temperature (< 10 9 K° ~ 100 KeV), neutron 
stars can be regarded as cold object. The main features of the theoretical 
models of neutron star matter will be summarized in Section 4. 

Given the EOS describig matter in the interior of a neutron star, its mass 
and radius can be obtained from the Tolman-Oppenheimer-Volkov (TOV) 
equations 



dP{r) 
dr 



= -G 



[e(r) + P{r)/c 2 } [M(r) + Anr 2 P{r)/c 2 
r 2 [1 - 2GM{r)/rc 2 ] 



(12) 



M(r) = 4tt / r' dr'e(r') , 
Jo 



(13) 



In the above equations, combining hydrostatic equilibrium with Einstein's 
equations of general relativity for a nonrotating star, P and e denote pressure 
and energy-density, respectively, and G is the gravitational constant. 

For any value of the central energy-density e c = e(r = 0) Eq.(|l2|) can be 
integrated outward until the star radius R, defined by the condition P(R) = 0, 
is reached. The star mass M(R) can then be obtained from Eq. Jl3]) . 
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Figure 4. Typical dependence of the neutron star mass (in units of the solar mass Mq) upon 
the energy-density at the center of the star (in units of the energy-density of symmetric 
nuclear matter at equilibrium). 

Fig.^| shows the typical behavior of the neutron star mass as a function 
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of the central energy-density. Denoting by e the central energy-density corre- 
sponding to the maximum mass, all stable configurations have e c < e c . 

The value of the maximum mass is mostly determined by the stiffness of 
the EOS, more incompressible neutron star matter (i.e. stiffer EOS) corre- 
sponding to larger M max . Threrefore, in principle measurements of neutron 
star masses may provide information on the behavior of the EOS at n > hq. 

The most precise experimental data—obtained from studies of the timing 
of radio pulsars, yield the average value □ 

M = 1.35 ±0.04 M , (14) 

thus constraining any realistic model of EOS to support a stable star with 
mass ~ 1.4 M Q . 

Neutron star masses can also be obtained from the analysis of binary 
systems containing an X-ray pulsar. A recent jdetermination of the mass of 
the Vela X-l pulsar yields M = 1.86 ± 0.33 M suggesting the possibility of 
a mass well above the canonical value of Eq.jl^). If confirmed, the existence 
of a neutron star with M > 1.8Mq would provide a very stringent constraint, 
ruling out theoretical models that predict very soft EOS. 

4 Theoretical models of neutron star matter at supranuclear 
density. 

Theoretical models of neutron star matter are based on different approaches, 
involving different degrees of freedom, which are expected to be applicable in 
different density regimes. 

In the range uq < n < 4hq matter in the interior of the star is believed 
to take the form of a cold uniform fluid, that can be described in terms of 
hadronic degrees of freedom, within the framework of either nonrelativistic 
many-body theorjfl or relativistic mean-field approaches!! 

At n just above tiq neutrons dominate, with a small fraction of protons 
and leptons in equilibrium with respect to the weak interaction processes (£ 
denotes either an electron or a muon) 

n^p + l + Vg , p + £—*n + ut. (15) 

For any baryon density n, minimization of the total energy-density with the 
constraints of baryon number conservation and charge neutrality fixes the 
relative aboundances of neutrons, protons and leptons. Note that, as the 
proton fraction is small, typically < 10%, theoretical calculations of the EOS 
involve the study of stronlgly asymmetric nuclear matter, whose energy is most 
sensitive to the isospin asymmetry term in Eq. (@) . The importance of isospin 
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asymmetry is illustrated in figj5| showing the calculated density dependence 
of the energy per particle of symmetric nuclear matter ((Z/A) — 1/2) and 
pure neutron matter ((Z/A) = 0). It clearly appears that unlike symmetric 
nuclear matter, whose energy exhibits a minimum, pure neutron matter is not 
self-bound at any density. 




12 3 4 

p/po 



Figure 5. Typical density dependence of the energy per nucleon of pure neutron matter 
(upper curves) and symmetric ppclear matter (lower curves), obtained from nonrelativistic 
nuclear many-body approaches^. The density is given in units of the equilibrium density of 
symmetric nuclear matter. 



As the density increases, different forms of matter, containing hadrons 
other than protons and neutrons, can become energetically favored. For ex- 
ample, the weak interaction process 

n + e — > £~ + v e , (16) 

leading to the appearance of a strange baryon, sets in as soon as the sum of the 
electron and neutron chemical potentials exceeds M s - (typically at n > 2n ). 
At larger density the production of A°'s is also energetically allowed. 

For any fixed baryon density, the relative aboundances of the different 
hadronic and leptonic species are dictated by the requirements of equilibrium 
with respect to weak interactions, conservation of baryon number and charge 
neutrality. Fig.^ shows theLparticle fractions resulting form the nonrelativistic 
many-body aproach of refJj. 
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Figure 6. Relative aboirmdances of hadrons and lcptons in a hadronic model of neutron star 
matter (taken from refja). 



The transition from nucleon matter to hadronic matter is associated with 
a softening of the EOS, as it basically amounts to replacing particles carrying 
large Fermi energies with more dilute, and therefore less energetic, strange 
baryons. Hence, using the EOS of hadronic matter as an input for the so- 
lution of the TOV equations leads to predict a value of the maximum mass 
significantly lower than that obtained using nucleon matter. 

At very large density (n > 4n ) a new transition, to a phase in which 
quarks are no longer clustered into nucleons or hadrons, is eventually expected 
to take place. In most theoretical studies of neutron stars structure, the EOS 
of quark matter has _heen estimated using the simple MIT bag model (for a 
recent review see refJiJ). 

The bag model assumes that quarks be confined to a finite region of space 
(the bag), whose volume is limited by a pressure B (the bag constant). Its 
simplest implementation includes only massless and noninteracting u and d 
quarks. Equilibrium with respect to the process (£ denotes either an electron 
or a muon) 



requires that the u and d chemical potentials satisfy (note that, as the mean 
free path of neutrinos in dense nucleon matter largely exceeds their typical 



d — > u + £ + vg , 



(17) 
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radius, neutron stars can be regarded as transparent to neutrinos) 

fid = Mtt + fit , fie = Hft, (18) 

whereas the constraint of charge neutrality implies 

2 1 

-n u = -n d + n e + n M , (19) 

where n q (q = u, d) and ng (I = e, fi) denote the quark and lepton densities, 
respectively. 

The energy-density reads 

E 3 



= (20) 



and the quark Fermi momenta are related to the total baryon density n 
through 

PF q = (* 2 f q n)^ , (21) 

j q being the number of quarks of flavor q per baryon. Finally, the pressure 
can be readily obtained from 

^=-(I)=-*+ii>^ < 22 > 

More realistic models include a massive strange quark and allow for one- 
gluon exchange interactions between quarks of the same flavor. At fixed 
baryon density the quark and lepton fractional densities are determined by 
the requirements of chemical equilibrium, charge neutrality and baryon num- 
ber conservation. Fig.|?] shows the results obtained including massless u and d 
quarks and strange quarks of mass m s = 150 MeV, neglecting all interactions 
and setting the value of the bag constant to B = 208 MeV/fm 3 . 

In the n — ► oo limit, the energy per baryon of nucleon matter, ejvM, grows 
linearly with baryon density, whereas in quark matter cqm ~ n 1 ^ . Hence, a 
transition from nucleon to quark matter is expected to take place in the inner 
region of neutron stars, provided the central density is large enough. 

In the simplest implementation, the phase transition is assumed to take 
place at constant pressure and chemical potential, so that matter density 
exhibits a discontinuity at the boundary between the two phases. According 
to this picture the star consists of a core of quark matter surrounded by 
nucleon matter. j—. 

The calculations of refJca, based on nonreativistic many-body theory for 
the description of nucleon matter and the bag model for the description of 
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Figure 7. Relative aboundanccs of quarks and leptons in /3-stable matter consisting of 
massless u and d quarks and strange quarks with mass m s = 150 MeV. All interactions 
between quarks are neglected and the bag constant is set to B = 208 MeV/fm 3 . 



quark matter, predict a transition region 5.4 < (n/no) < 9.8 for B = 200 
MeV/fm 3 and 4.9 < (n/no) < 7.5 for B = 122 MeV/fm 3 . 

The transition to quark matter makes the EOS softer, thus lowering the 
value of the maximum neutron star mass. The authors of ref.Eil report a 
reduction AM max of - 10 % and - 20 % for B = 200 and 122 MeV/fm 3 , 
respectively. 

5 Will detection of gravitational provide new insight ? 

As the pattern of nonradial oscillations of neutron stars depends upon the 
EOS describing matter in the interior of the star, detection of gravitational 
wave emission associated with the excitation of these modes can in principle 
provide additional constraints on theoretical models of the EOS. 

For example, the complex oscillation frequencies of the axial (i.e. odd par- 
ity) modes of a nonrotating star are eigenvalues of a Schrodinger-like equation 
whose potential Vt(r) explicitely depend upon the EOS according toE3 

e M»0 

V e (r) = {1(1 + l)r + r 3 [e(r) - P{r)} - 6M(r)} , (23) 
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where 

*L = I ^ (24) 

dr [e(r) + P{r)\ dr V ; 

Using realistic EOS one obtains strongly damped eigenmodes, called w-modes, 
whose frequencies exhibit the behavior diaplayed in figM 




Figure 8. Frequencies of the first axial w-mode of a nonrotating neutron star, plotted as 
a function of the star compactpass (M/R). The different curves correspond to different 
models of EOS (taken from ref.F 2 !). 



The pattern of frequencies strictly reflects the stiffness of the different 
EOS in the relevant density region (typically 1 < (n/no) < 5, softer EOS cor- 
reponding to higher frequencies. For example, the curve labelled B has been 
obtained from a model including nucleons, nucleon excitations and strange 
baryons, predicting a very soft EOS and a maximum mass of f .42 M®. 

It is interesting to note that the dependence upon the the ratio (M/R) is 
rather weak, so that detection of a given frequency may allow one to discrim- 
inate between different models of EOS (e.g. between model B, corresponding 
to hadronic matter and model WFF, correponding to nucleon matter) regard- 
less of the compactness of the star. 

A different class of nonradial oscillations, the g-modes, are associated with 
the occurrence of a density discontinuity in the interior of the star. Early stud- 
ies of these modes focused on discontinuities that are known to be produced 
by the changes of chemical composition in the low density region of the crust, 
corresponding to a fractional distance from the surface (d/R) < 10 %. 
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As pointed out in Section 4, a different discontinuity may be produced by 
the transition from nucleon to quark matter, which is expected to take place 
at much larger density and involve a much larger density jump. 

It has to be emphasized that, as nucleon and quark matter need not to be 
considered independently charge neutral, in general the deconfinement tran- 
sition may involve the occurrence of a mixed phase (e.g. bubbles of quark 
matter in nucleon matter at lower density or viceversa at larger density) ex- 
tending over a sizable region of space, rather than a sharp separation between 
the two phasesEl. Whether the mixed phase is energetically favored depends 
upon the balance between the gaixLin volume energy and the loss associated 
with Coulomb and surface energyc^l. If Coulomb and surface energies prevail, 
a density discontinuity is expected to occur. 
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Figure 9. Frequencies of the f- and g-mode of a neutron star of mass M = 1.4Mq as a 
function of the average density. Solid, dashed and dot-dash lines correspond to density 
jumps of 10 %, 20 i%iand 30 %, respectively. All stellar models have polytropic exponent 
T = 2 (see Eq.@)EJ 

The authors of refill have recently carried out an exploratory calcula- 
tion of the frequencies of both the fudamental f-mode and the g-modc of a 
nonrotating neutron star, using a simple polytropic EOS with a density dis- 
continuity An located at density n = no- 

[(l + An) T n T n <n D 
Kn n > nrj + An . 
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The results of refJia are summarized in fig]9|, showing the f- and g-mode 
frequencies as a function of the average density of the star, whose mass is 
kept fixed at the canonical value of 1.4 M©. It clearly appears that, while 
the f-mode frequency exhibits a nearly linear growth, largely unaffected by the 
abrupt structural change associated with the phase transition, the frequencies 
of the g-mode strongly depend upon the jump in density An. Hence, a simul- 
tanuoeus measurements of the frequencies of the two modes would provide 
information on both size and location of the discontinuity. 



6 Conclusions 

The extent to which the currently available empirical information constrains 
the EOS of cold /3-stable strongly interacting matter is illustrated in fig JlOj. 
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Figure 10. EOS of cold /3-stable strongly interacting matter resulting from a variety of 
different theoretical approaches, based on both hadronic and quark degrees of freedom. The 
shaded areas show the Feeion forbidden by causality and that corresponding to M m ax < 
1.4 Mq (taken from ref£j)- 
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Although all EOS resulting from nonrelativistic approaches are known to 
be plagued by a noncausal behaviour at very large density, it appears that the 
constraint P > e is fulfilled by all EOS over a wide range of energy-density. 

The requirement that the maximum mass be larger than 1.4 Mq does not 
provide a stringent constraint either, as all reasonable models support a stable 
star with M — 1.4 Mq. However, as pointed out in Section 3, the situation 
may change should the mass of the Vela X-l pulsar be confirmed to be larger 
than 1.8 Mq. In this case the shaded region in the lower right corner of fig.|lO| 
would considerably extend upward, thus ruling out models that predict very 
soft EOS. 

The observation of gravitational waves emitted by neutron stars may pro- 
vide valuable new insight, allowing to further constrain theoretical models. 
Detection of gravitational radiation in the relevant few KHz frequency range 
may become possible with the recently proposed gravitational laser interfero- 
metric detector EURO E3. 

As far as theory is concerned, it has to be mentioned that the recent 
progress in Quantum Monte Carlo has made it possible to apply this approach 
to the study of the EOS of pure neutron matter E3. The extension to the case 
of /3-stable matter, that appears to be feasible, would be of great importance 
for the understanding of neutron star matter in the region of not too high 
density (n < 2n ), where the use of nucleonic degrees of freedom is likely to 
be reasonable. 

At larger densities a transition to quark matter is expected to take place. 
A better theoretical description of this region will require a careful analysis 
of the nature of the phase transition as well as a more realistic description of 
the deconfined phase. In this context, a pivotal role is likely to be played by 
detailed quantitative studies of the possible occurrence of the color supercon- 
ducting phase predicted by the fundamental theory of strong interactional! 



References 



K. Huang Statistical Mechanics (Wiley, New York, 1963). 

H. Heiselberg, talk given at the Workshop on "Compact Stars in the QCD 



Phase Diagram", Copenhagen, August 15-18, 2001. |astro-ph/0201465 
S. Shlomo and D.H. Youngblood, Phys. Rev. C 47, 529 (1993). 
Ya. B. Zcl'dovich, Sov. Phys. JETP 14, 1143 (1962). 
S.E. Thorsett and D Chakrabarty, Ap J 512, 288 (1999). 
M.H. van Kerkwijk, |astro-ph/0110336 . 

H. Heiselberg and V.R. Pandharipande, Ann. Rev. Nucl. Part. ScL, 50, 
481 (2000). 



eos: submitted to World Scientific on February 8, 2008 15 



8. N.K. Glendenning, Compact Stars (Springer, New York, 1997). 

9. I. Vidana et at, Phys. Rev. C 62, 035801 (2000). 

10. H. Heiselberg and M. Hjorth- Jensen, Phys. Rep. 328, 237 (2000). 

11. A. Akmal, V.R. Pandharipande and D.G. Ravenhall, Phys. Rev. C 58, 
1804 (1998). 

12. O. Benhar, E. Berti and V. Ferrari, Mon. Not. R. Astron. Soc. 310, 
797 (1999). 

13. H. Heiselberg, C.J. Pethik and E.F. Staubo, Phys. Rev. Lett. 70, 1355 
(1993). 

14. G. Miniutti, J. A. Pons, E. Berti, L. Gualtieri and V. Ferrari, paper in 
preparation. 

15. F. Weber, Pulsars as Atrophysical Laboratories for Nuclear and Particle 
Physics (IoP Publishing, Bristol, 1999). 

16. EURO: Europe's T hird Generation Gravitational Wave Observatory. 
Proposal available at http://carina.astro.cf.ac.uk/geo/euro J (May, 2000). 

17. K.E. Schmidt and S. Fantoni, Phys. Lett. B446, 99 (1999). 

18. M.G. Alford, K. Rajagopal, S. Reddy and F. Wilczek, Phys. Rev. D 64, 
074017 (2000). 



eos: submitted to World Scientific on February 8, 2008 16 



